Diagram expansions in classical stochastic field theory 



L. I. Plimak, M. Fleischhauer* , M. J. CoUett and D. F. Walls 
Dept. of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand 
- - - * Sektion Physik, Ludwig-Maximilians Universitdt Miinchen, D-80333 Miinchen, Germany 

' (February 7, 2008) 

0\ ; 

On , A diagram approach to classical nonlinear stochastic field theory is introduced. This approach 

is intended to serve as a link between quantum and classical field theories, resulting in an indepen- 
^ ' dent constructive characterisation of the measure in Feynman path integrals in terms of stochastic 

5^ , differential equations for the paths. 
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I. INTRODUCTION 



o 

^ , Diagrams are commonly associated with either Feynman's work on quantum electrodynamics, or with Matsubara 
and Keldysh techniques in statistical mechanics, so that diagrams appear to be something to do with subtle problems 
in relativistic quantum field theory or the theory of phase transitions. This view is totally misleading. The concept 
^ . of dia grams is very general and can be applied to virtually any problem in physics. By encoding details of a physical 
c/2 problem as expressions for the lines and vertices, a solution to this problem is readily produced following general 
diagram rules. 

There exists an additional dimension to the diagram approach not yet fully appreciated, namely its ability to 
y . establish links between approaches based on inconsistent techniques, like quantum and classical field theories (QFT 
K • and CFT). The diagram approach to the QFT is common knowledge. In this paper, we develop a diagram approach 
^ ' to the theory of stochastic differential equations (SDEs), which is the formal ground of the CFT. The central formal 
Q result of the paper is a general relation which we call the causal Wick's theorem. It plays the same role in deriving 
O diagram series for SDEs as Wick's theorem in QFT plays in deriving the Feynman (say) diagram series. We then 
find that any SDE is formally solved by diagram series of a certain structure, which we call causal diagram series. In 
I ' these series, the propagator is determined by the linear part of the SDE, and vertices by the nonlinearity and noise 
^ , sources. This result may be interesting in itself since it allows one to extend methods specific to diagram approaches 
(Nj ' (e.g., the Dyson equations) to SDEs. 

On ' The formal relation between diagram series in the QFT and in the CFT will be investigated in papers to follow. 
Here we note only that the (pseudo-)stochastic measure on the solutions of an SDE is exactly the measure in certain 
Feynman path integrals: solutions to the SDE are the Feynman paths. This may result in better understanding of 

the path-integral approach and also in the possibility of characterising path integrals constructively, in well-defined 

' mathematical terms, and, indeed in a way independent of perturbation theory. 

, From an even more fundamental perspective, we encounter a very important link between regularisations in the 
' diagrams and stochastic calculus in SDEs. Technically, we find that regularisations in a diagram series which make 
^ , diagrams convergent affect the respective SDE in such a way that it becomes mathematically defined within the normal 
I ' calculus. We also find some indication that a renormalisation procedure, which formally consists of (i) regularisation 
' O ' and (ii) a limiting procedure removing this regularisation in a certain way, appears from the point of view of the 
^ •■ respective SDE as a limiting procedure specifying the stochastic calculus. This link between renormalisation and 
stochastic calculus suggests more extensive mathematical investigation, which may lead to better understanding of 
the formal grounds of quantum field theory. 



II. CLASSICAL STOCHASTIC SELF-ACTION PROBLEM 

A. Stochastic calculus and regularisations in SDEs 

We consider a c- number field ip{r, t), which satisfies a generic equation with a source s(r, t), 

C^{r,t)^s{r,t), (1) 

where £ is a differential operator. The source s is dependent on the field ^ and this dependence may include 
randomness; that is, in general (|l|) is a stochastic diflFerential equation for the field t)- To be specific, we assumed 
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that -0 is a classical field in one-dimensional space; to consider other situations, e.g., a multi-mode field, or a field in 
three-dimensional space, one should simply replace in the relations below J dr by J2r (^■^■' then a mode index), 
or / (fir, respectively. For a single-mode field, the variable r and the summation should simply be dropped. We shall 
often resort to this last case for simplicity of examples. 

Formally solving (|l|) turns it into an integral equation, (assuming s(r, t) — > as t — > — oo) 

i'ir, t)= I dt' [ dr'Gir, r',t- t')s{r', t') + Mr, t), (2) 



where G(r, r' ,t — t') is the retarded Green's function of the equation (|l]), 

CG{r,r\t-t') = &{r-r')5{t~t'), G{r,r' ,t - t') = Q,t < t' , (3) 

(so that integration in (j^) is in fact from minus infinity to t), and V'o is the field, i.e. iJj — ipo before the source is 
on. It obeys the free version of the equation (0), 

CMr,t)^0. (4) 

Note that the existence of a retarded Green's function is not guaranteed for an arbitrary £; conditions (^ single out 
equations that may have physical meaning. 

Although these definitions (as well as results below) are quite general, three cases are of major practical importance: 
the nonrelativistic Schrodinger equation, 

C^idt + \dl, (5) 

the wave equation, 



L = di-dt, (6) 



and the relativistic Klein-Gordon equation. 



c = d^- dl + 1, (7) 

where units were chosen so as to remove dimensional constants. Spatial Fourier transformation, 

i,(r,t) = J dke^'^'^Mt), (8) 

reduces these equations to their single-mode cases: 

C = tdt-'^, Gk{t) = -ie{t)e-^^'''\ (9) 

C = d^ + k^, Gk{t) = e{t)fimkt, (10) 

C = df + k^ + l, Gfc(t) = 6'(i)sintA/FTT, (11) 

where the corresponding Green's functions are also given. 

These Green's functions all have a singularity at zero time, so problems arise if the source s is a singular function 
as well. E.g., ii C — idt — fc^/2 and s contains white noise, definingmathematically the inhomogeneous Schrodinger 
equations (|l|) and (|^) requires specification of a stochastic calculus |]l| . Another way around this problem is found by 
assuming that the Green's function and/or the source are regularised, e.g., properly smoothed, and a corresponding 
limiting procedure is employed for restoring their singular values. The choice of stochastic calculus then becomes a 
result of this regularisation. For example, regularising the source by introducing finite correlation times means that 
normal calculus holds at all stages of the regularisation procedure, so that in the end the Stratonovich calculus is 
recovered. 

As we shall argue now, regularising G may lead to Ito calculus. Consider, for example, the equation, 

tdtm = emw'{t), (12) 

where e is a constant and W' is the derivative of the Wiener process. It can only be defined in the sense of the theory 
of generalised functions, 
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j dtW'it)ip{t) = - j dtW{t)ip'{t), (13) 

where ip{t) is a "good" function. It would suffice for to have finite support and be continuously differentiable, so 
that <f' is a function of finite variation. It is then clear that for ipW' to be defined, it would suffice for -0 to be 
continuously differentiable, but this is certainly not consistent with Eq. (p^). 

The standard way around this problem is to define / ipW'dt — J 'ipdW as a stochastic integral so that the equation 

ij) = — ie TpdW is defined. Instead, we replace an undefined integral equation, tp = tpo — ie tpW'dt', by a 
regularised integral equation, 

i^it) = + / Greg(t - t')^{t')eit')W'{t')dt'. (14) 



Grcg(t) is a regularised retarded Green's function, which is (i) causal, Grcg{t) = 0,t < 0, (ii) a given number of times 
continuously differentiable, and (iii) in a certain sense close to G{t). For example, {k > 1) 

G,cg(t)--ie(t)(l-e-")'+' (15) 

is k times continuously differentiable. F here is an arbitrary positive parameter, and the final limit F ^ cxd is implied. 
We have also introduced a truncating factor e{t) into the white noise factor W': e{t) < £ is infinitely differentiable, 
and a negative T' exists such that e{t) = e if t > T' and e{t) = if t < 2T'; this is necessary to have a consistent 
m-field formulation. Eq. (|l^) is now consistent with the assumption of ip being k times continuously differentiable. If 
this is indeed the case, the factor Gicg{t — t')%p{t')e{t') at a given t > 2T' is (at least) continuously differentiable and 
has a finite support {t' : 2T' <t'< t}, so that the integral on the right of (Hj) is defined; then, continuity of its A:-th 
derivative by t follows from the fact that W' is, loosely speaking, no more singular than a 5-function. (For t < 2T', 
( p^ ) reduces to ip = ipo-) 

An indication that in the limit F — )■ cxd a solution to Eq. (|lj) approaches a solution to the Ito differential equation 
( |l^ ) may be seen from the following considerations. As a generalised function, e{t')W'{t') may be approximated by 
a discrete sum of J-functions, 



k=l 

where tk = 2T' + fcAt, and Ai is a discretisation scale. Eq. ( p^ then has a unique solution, 

oo 

=i^O+Y. ^kGregit - tu)^{tk), (17) 
fc=l 

where ipltk) may be found recurrently, 

m — 1 

^{tm) = -00 + X! ^kGrcgiUn - tk)lp{tk). (18) 



k=l 



The sum in (|7|) is in fact finite, so that 0(t) inherits all the "goodness" of Grcg(i)- It is now easy to see that if 
FAt ^ 1, the integral in ( |l4[ ) coincides with the partial sum —i^'ip{tk)s{tk) \W{tk+i) — VF(tfc)]; so that in the limit 
At ^ we recover the Ito integral —i J^j,, dW{t')iJ;{t')£{t'). From the practical point of view, this "proves" the 
hypothesis that in the limit F — » oo the Ito calculus is recovered, since any real calculation implies time discretisation. 
In order to prove it mathematically, one should commute the limits: first J2T=i ^kS{t — tk) e{t')W'{t') and second 
F — > oo, while the above considerations imply the opposite order of the limits. 

Below we shall see that G(r, r' ,t — t') becomes the propagator in the diagram series so that regularising it prevents 
ultraviolet divergences in the diagrams (this is nothing but the Pauli-Villars regularisation known in QFT [||). The 
choice of stochastic calculus is hence connected to regularisations in the diagram series and hence to renormalisations. 
Investigation of this connection in full is a major undertaking; it should also include the continuous-space limit which 
we do not even try to tackle. At this stage, we just assume that G and/or s are regular (or regularised if necessary), 
so that all our expressions make sense. Note that this means that the integral equation (^) rather than the differential 
equation (Q) is considered: the latter only emerges when regularisations are removed. 
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B. Characterisation of the stochastic sources and fields 



In order to have a fully defined stochastic differential equation, the functional dependence of the source on the full 
{local or microscopic) field, t), is needed. Such a charcterisation of the source is usually given by the physical 
model considered. Then to solve equations (|l|) or (0) means finding the field -tp as function of the m-field ipoij-jt). 

For example, in a homogeneous linear medium, s{r, t) = X'0(^: where x is the linear susceptibility of the medium. 
In general, when randomness is also included, the source is a stochastic variable whose properties are determined 
by a probability distribution, P(s[?/'), of the function s conditioned on the full field ip(r^t). It will be convenient 
to distingish between real and complex fields. We will start from the former case and will generalize the results to 
complex fields at the end. We introduce a characteristic functional S{a\ip) corresponding to P(s|V') as, 



S{a\i;) = e'-% = J e"^P(s|i^), 



(19) 



Here, a{r,t) is an arbitrary "good" real function and / Ds denotes a functional integration over the "trajectories" 
s{r,t). We use an abbreviated notation in which as denotes J dxa{x)s{x), where x = {r,t} and J dx = J drdt. 
Differentiating S over a{x) produces multi-space-time averages of the source, conditioned on the full field, 



s{xi) ■ ■ ■ s{Xn) 



Qn 



ip da{xi) ■ ■ ■ da{xn) 



S{a\ilj) 



(20) 



To specify the functional dependence of the source on ij: it is convenient to introduce generalised susceptibilities (or 
just susceptibilities) x'"™'"'' (^i' ' ' ' j ^^m; a;']^, ■ • ■ , a;^). They are coefficients in the series expressing cumulants of the 
local source s{x) in terms of the powers of the local field ipi^)- 



{x) =x''^°^W+ dx'x^^^^\x-x')4>{x') 



+ i I dx'dx"x^^'^Hx;x',x")^{x')i^{x") + 



s{x)s{x') = s{x) s{x') +x^^'°\x,x')+ I dx"x^^'^'>{x,x';x")i:{x") + 

1p 'ip tp 



(21a) 
(21b) 



With them the generating functional can be given a simple form 

oo 

5(a|V')=exp V a™x^"'"V", (22) 

^-^ mini 

m,n— 

where 

dXfjidx ' ' ' dXj^ X 

a(x„)x(™'") (xi, • • ■ , x'l, • ■ • , x'Mx[) ■ ■ ■ ^«). (23) 

The susceptibilities should obey the causality condition, 

X^'"'"Hxi,---,a;„;x;,---,0 =0, max{t[, ■ ■ ■ ,0 > max(ii , • • • , t„0 ; (24) 

this condition can also be formulated as the latest argument of a susceptibility is always an output one; hence x^°'"^ = 
for all n. 

The quantity x'-^''^-'(a^) is a given source (non-random). For m = 1 we find susceptibilities proper: the linear one, 
■0'-'^'^^ and the nonlinear ones, ^/j^^'"' for n > 1. If x'™'"-* = for m > 1, the dependence of the source on the local field 
is not stochastic, and equation (|l]) is not stochastic either; it is linear, if x*-^'"'' = for n > 1, and otherwise nonlinear. 
Non-zero x^™'"^ for m > 1 introduce stochasticity: non-zero x^"^'^^ are cumulants of a given random source, while 
^(m,n) both m > 1 and n > describe how the statistics of the source depends on the field. 

In practice, as a rule, only a finite number of generalised susceptibilities are non-zero. To be specific, we restrict 
ourselves to the case of only x^^'^\ x'"^'^\ x'^'^\ ;^(2'0) g^jj(j x'^'^^ non-zero; these are exactly the susceptibilities shown 



a™x'™'"V" 



dxi 
a{xi) ■ 
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explicitly in Eqs. (^l]). This case corresponds to a source which is Gaussian if conditioned on the full field (this 
reservation is important). Also, the susceptibilities are commonly local: e.g., x^^'^H^!^') ~ ^i^ ~ so that the 
source at the point x depends only on the local field at the same point. The causality condition ( |24| ) is then satisfied 
automatically. Locality gives more physical sense to the equations, but leads to mathematical problems. In this paper, 
we assume that these problems are overcome by regularisations. 

Solving equations (jl) and (|^) means finding the field dependence on the in- field V'o- It is immediately obvious that 
this is equivalent to finding the dependence of the source on the in-field, i.e. given P{s\^) finding the probability 
distribution n(s|^/'o). Then, once the "macroscopic" probability distribution n(s|'(/'o) is known, solving for the field 
is straightforward. To characterise the field, we introduce the characteristic functional of the multi-space-time field 
averages. 



$(C) = $(C|V'o) = 



•00 



so that 



•(/'(xi) • • ■ ■0(a;„) = tpixi) ■ ■ ■ ijj{x„) 



gn 



(25) 



(26) 



Vo dC{xi) ■ ■ ■ dC{xn) 
Since V' — Gs + ipo, the field statistics are effectively that of the source, and we have, 

$(() = J L>se'^(^"+'^°)n(s|Vo) = c';*'E(CG|V'o), (27) 

where [CG'](a;) — J dx' C,{x')G{x' ,x). Here we have introduced the characteristic functional S(a|V'o) corresponding to 
n(s|Vo): 



E(alV'o) = = j Dse^'n{s\i:o). 



(28) 



In the complex case, our definitions must be generalised to allow one to consider averages containing the fields and 
sources as well as their complex conjugates. Thus the notation 11(317/10) used for the real case, must be replaced for 
complex fields by n(s, s* j^/'o, "^o) 5 the functional integration in the complex case will be denoted as J D^Dip* and 
J DsDs*, etc. The characteristic functionals S, E and $ are defined as. 



S'(a,aT|i/',V*) = / DsDs*c' 



(29) 
(30) 
(31) 

where q;(x), (x), C(a;) and C'^ix) are arbitrary "good" functions (^ here is just a notation, not relevant to Hermitian 
conjugation). One can also consider a, and (, &s pairs of complex-conjugated functions, at = a* , = C* ■ Then, 
e-g-, 



S(a,«t|^„,^*) = J i?.si?s*c"'^+"^'n(s,s*|Vo,V'o): 



S*{xi) ■ ■ ■ S*{Xm)s{x[) ■ ■ ■ S{x'^) 



d 



im+n 



da{xi) ■ ■ ■ da{xm)da^ (x'-^) ■ ■ ■ da^{x'j^) 
xS{a,a^\i},'il)*)\a=a^=o. 



(32) 



The rest of the above relations change accordingly. 



C. Relation between P(s|'!/)) and n(s|'!/)o) 

Our goal now is to find a formal solution to Eq. (^ in the form of a relation between the "microscopic" and 
"macroscopic" probability distributions, P(s|'0) and n(s|-!/)o)- Assume discretisation of the time axis, s{t),ip{t) 
s(ife), V'(^fe)j where = kAt, k = —00, • • • , 00; a final limit of A< — > is implied. We omit the spatial variable as 
irrelevant. With discretisation, averages of the source at the given field are given by a functional integration {t, t' ,■ ■ ■ t" 
are among tk), 
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/oo 
sit)s{t')---sit")Pis\ij) rr dsitk 



(33) 



k— — oo 



The distribution P(s|'0) is causal, ie, s(tfe) depends on il^itm) only for m < k (which is equivalent to the above causality 
condition for susceptibilities). This allows one to introduce reduced probability distributions, 



-Pm (S<t„ |V'<t„) 



/oo 
p(siv) n ^^(^'^) 
; — ™ I 1 



k—m+1 



where 



*<t,r. = {•s(^fe) : fc < m} , 
V'<t,„ = {V'(^fc) : k <m} . 



(34) 



(35) 
(36) 



The fact that Pm depends on V'<t,„ and not on the whole i/j leads to the following causality condition for the averages, 

d 



^in ^ s{t)s{f)---s{t") =0, t„ >max(i,t',---t"). 

Oip[tm) 4> 

The unravelling of the source statistics in time is described by the conditional probability distribution, 

Pm (s<t„ |V'<t„) 



Pm (s(tm) |s<t„_i,l/'<t™) = 



Pm-1 {s<tm-i |V'<t„-i) 



(37) 



(38) 



In turn, 



Prn(s<tr^\tp<t„J = Y]. Pk {s{tk)\s<t^_i,1p<t,,) 



(39) 



k— — OD 



Since we always observe a system only for finite times, P,,,. with m large enough is actually as good as P(s|V'), so that 
we can write 



p(siv)= n Pk{sitk)\s<t,_,,i^<t,] 



(40) 



fc— — oo 



(More rigorously, this relation implies w-limits for both plus and minus infinity.) 

The advantage of the "unravelling" representation (p]| ) for P{s\^) is that it is perfectly designed so as to accept 
the dynamical relation between the source and the field, Eq. (^. With time discretisation, it is understood as. 



(41) 



Note that since G{t) is both causal and regular (or regularised), G{0) ~ and hence the latest source value to 
contribute to ipltm) is s(i„i-i), ie, i^{tm) depends on s<tm-i (V'o is either non-random or uncorrelated with s). Then, 
both s<tm-i and V'<i,„ in the conditional probability P™ (s(im) |s<t„_i , V''<tm) are effectively non-random and we 
can write 



n„i s{t 



s<t„ 



-i>(V'o)<t^) = Pm (s(i„) s<t^_i,(Gs-|-?/;o)<t^) 



(42) 



Here, !!„ {^s{tm) s<t,„_i, (■0o)<t ^ is the probability distribution for the source at t™, conditioned on its own prehis- 
tory and the in-field. Relation (B3) solves the self-action problem for the source, expressing its actual (macroscopically 
observable) statistics in terms of the (microscopic) relations characterising the system. For the multi-time probability 
distribution for the source, conditional on the in-field, we have. 



n(s|V'o) = W n„ (s(t™) s<t,„_i, (V'o) 



(43) 



rn— — OO 
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and hence, 

n(s|Vo) = P(s|Gs + V^o). (44) 

This relation originates in Eq. (^) (regularised!) and causality. 

The above derivation does not depend on whether we deal with a real or complex field. However, following the 
above convention regarding notation, in the complex case it will be written as 

n(s, s* IVo, V'o) = Pi^. ^IGs + V^o, G*s* + ro)- (45) 
Relation (Q) applies in the real case. 



D. Causal Wick's theorems 

Our aim is now to rewrite relation in terms of the characteristic functionals. To this end, we first rewrite this 
relation using the shift operator, 

n(slV'o) = exp (^Gs^ P(s|^)|v,=^o. (46) 

Here, we again use condensed notation, -§^Gs = J dxdx' -q:^^G{x;x')s{x'); note that the spatial dependence of the 
fields is restored. Then, 

E(a|i/'o) = J Dsexp ^as + '^^^^ ^(s|V')lv='/'o 

= "''P (^^^) / Dse'^^PisM^^^o = exp (^g|-) 5(a|V)lv>=Vo- (47) 

It is worth noting why we have to write, ^(' ■ ■)\jp=ipo, rather than just The problem is that V'o is a solution 
to a free equation, ie, is a derivative with constraints, whereas when applying the shift operator, exp (^-^Gs^ , to 

P{s\ip) in order to turn it into P{s\Gs + V'), one has to assume that i/; is arbitrary. Since P{s\tp) is indeed defined for 
an arbitrtar y ih , the whole situation is consistent. 

Relation ( p7| ) is startlingly reminiscent of relations known in QFT that express Wick's theorem for bosonic operators 
as a differential operation ||^,^ , and we shall call it the causal Wick's theorem. Below we shall see that it plays the same 
role in deriving diagram series for the SDEs as Wick's theorem proper plays in deriving diagram series for interacting 
bosonic fields. The assumptions it is based on are: (i) causality, (ii) Eq. and (iii) the condition G(r, r',0) ~ 0. 
Note that the last condition is meaningful only for regularised SDEs. 

For the functional $(C), relations ( |47| ) and ( p7|) result in, 

*(C) = exp (^G^) exp (C^) 5(a|V')|a=o,v=^o- (48) 



dt/j da 

This relation will be referred to as a generating formula for causal diagram series. 

In the complex case, relations expressing E and $ over S (i.e., the causal Wick's theorem and the generating formula 
for causal series) are found to be, 

/" d d c) \ 
S(a,at|V'o,^o) = exp (^—G-^ + ^G* —j ^(a,a^|i^, Vt)U=^o,v,t=^S, (49) 

In these relations, ip^ V'^ is a pair of "good" functions which can be either arbitrary or complex conjugated. We have 
assumed that 5(0;, IV', V'*) may be regarded as an analytic function of separately '0 and '0*; in all practical examples 
S = exp(Polynomial of ^, '0*) (cf Eq. pTl), so that this assumption is valid. Note, without going in detail, that 



assuming V'^ "good" is consistent only with a regular G{x]x'), so that relations (49) and ( |50| ) imply regularisations; 
this certainly applies to (E^) and (El) as well. 
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III. DIAGRAM EXPANSIONS AND SDES 



In QFT, relations similar to the causal Wick's theorems are used in order to express Wick's theorem proper in 
a compact form as a differential operation Then, if this form of Wick's theorem is plugged into the standard 

perturbation approach relations for Green's functions of the interacting bosonic fields are found which are 

structurally identical to relations ( ^8| ) or ( |50| ) for classical averages (this identity is straightforward for non-relativistic 
quantum fields [^[0], but becomes more involved in relativistic QFT j^). In principle, this identity is all that is 
necessary in order to establish the link between the quantum and classical field theories. For this purpose, diagram 
series (which both in QFT and CFT may be found by expanding generating formulae in power series) are redundant. 
Moreover, establishing the quantum/classical link in terms of the diagram series, strictly speaking, degrades the 
approach since it brings over the problem of convergence of the diagram series and confines the link to the perturbation 
theory (whereas, e.g., for the anharmonic oscillator linking the generating formulae leads to exact results |^,||). 

Despite all this, we believe introducing diagram series in CFT may be beneficial. First, diagrams are common 
language in QFT, and introducing them in CFT eliminates the seeming incompatibility between the q- and c- number 
techniques. Second, diagram expansions in CFT are not introduced as a computational tool. They are regarded 
only as a leading consideration, visualising certain structural properties of an SDE and/or the respective generating 
relation. These properties should then be proved directly. (Similarly, in QFT analysing diagram series leads to 
relations for observable quantities, like Dyson equations. As a rule, these relations may be derived independently.) 
Hence rigor of the diagram approach is not an issue. Third, in view of the encountered link between regularisations 
and stochastic calculus, introducing diagrams related to SDEs becomes a must: after all, those are diagrams that 
diverge! 



A. Causal diagram series 



For simplicity, we confine ourselves to a real SDE, with the noise source specified by the set of susceptibilities shown 
explicitly in Eqs. (pT|). Examples of series related to complex SDEs will be presented elsewhere Q. 
Graphically, the quantities ipoix) and G{x;x') will be denoted as lines. 



G{x;x') = ^^-—^'j, (51a) 
Mx) = I — 1 , (51b) 



while C(a;) and the susceptibilities as vertices, 

C{x) = I " • } , (52a) 



X(i^0)(:,) = |._^^|, (52b) 
X^^'^\x;x') = l^ "'"^ " |, (52c) 

X^'''^ix;x',x")^l\'\ "I, (52d) 



X(^'")(x,x')= , (52e) 



X^'-'\x,x';x")=^ j- (52f) 

The "time arrow" , which is drawn below each diagram, distinguishes graphically input and output arguments of 
the lines and vertices: The argument of tpoix) is regarded as a line output, as well as the "future" argument (i.e., 
x) in G{x]x'); the "past" argument (i.e., x') in G{x;x') is regarded as a line input. A generalised susceptibility. 
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X 



), is a quantity with n vertex inputs and m vertex outputs. By definition, the argument 
of C{x) is a vertex input. 

We now expand all exponents in the generating relation ( ^ ) in power series, and consider a certain term in these 
series. (In detail, this procedure is discussed in the appendix.) We see that a{x) always occurs convolved with a 
vertex output, while g^^jy is always convolved with a line input. Since q2{x') ~ S{x — x'), differentiating by a's leaves 
all vertex outputs pairwise convolved with line inputs. No free vertex outputs or line inputs may remain; terms with 
unequal number of these give zero. Similarly, derivatives g^jy leave vertex inputs convolved with the line outputs of 
the propagators G(x; x'); "surviving" ip^s become tpa^s. As a result, we find all vertex inputs pairwise convolved with 
line outputs, and no free arguments remain. Graphically, convolved input-output pairs are denoted by connecting 
respective ends of the lines to the vertices, e.g.. 



(53a) 
(53b) 




dxC{x)ijjo{x), 

dxdx'C{x)G{x;x')x^^''^Hx'), 

= J dxi ■ ■ • dx5({xi)C{x2)G{xi; X3)G{x2; X4) x 

x'^'^'^\x3,X4;x5)il;o{x5), 

dxi ■ ■ ■ dx(iC,{xi)G{xi\X2)x^^''^\x2;x3,Xi) x 



G{x3;xz)G{xi;x(i)x'''^'°\x,i, x&). 



(53c) 



(53d) 



Expressions of such structure will be called causal diagrams. Note that the "time arrow" applies to all elements in a 
diagram. 

In general, a causal diagram is a product of the basic elements — lines and vertices, where some line inputs and 
outputs are pairwise convolved with, respectively, vertex outputs and inputs. A diagram containing free arguments is 
called incomplete] e.g., (^ ) and ( ^2|) are legitimate incomplete diagrams. A diagram without free arguments (cf (|5^)) 
is complete. As is shown in the app endix, th e functional <&(C) is expressed as a sum of all complete causal diagrams 
which may be built of the elements (51a -521), with certain coefficients; the rules for calculating these coefficients may 
be found in textbooks. Incomplete diagrams occur, e.g., in diagram expansions for the field averages. Formally, thes e 
are found by "stripping" the C vertices from the complete diagrams. The simplest example of such a diagram is (|5Hj): 



it contributes to i}}{r, t) and is found by stripping the C vertex from (53a). 

It is also worth noting that whereas the "time flow" in diagram s is from left to right, in analytical expresssions, 
as a rule, time increases from right to left. Eg, in the diagram in (53b) the C, vertex is on the right, whereas in the 
analytical expression ^Gx'^'°^ — J dxdx'({x)G{x;x')x^^'^Hx') the natural position of the function ((x) is on the left. 
The diagram notation thus commonly has an inverted order of objects compared with the analytical notation. 



B. Connected diagrams and field cumulants 



A diagram that graphically consists of a number of separate sub-diagrams without common elements is called 
disconnected; otherwise, a diagram is connected. A disconnected diagram is a product of its connected components. 
In the above examples, all diagrams are connected; however, full causal diagram series for the functional $(C) contain 
all possible complete diagrams, connected as well as disconnected. 

To get rid of disconnected diagrams, one should describe the field in terms of its cumulants rather than averages. 
Formally, the field cumulants G^'^\ n — 0,1,- ■ ■, are defined as, 

oc 

$(C)-exp^-C"C("), (54) 

n—l 

where (n = 1, 2, • • •) 
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C"C(") = J dxidX2 ■ ■ ■ dXnC^"\xi,X2, - ■ ■ ,Xn)C{xi)aX2) ■ ■ ■ C{Xn). (55) 

In particular, 

^4J{xj =C^^\x), (56) 

^(x)V'(x') = C(2) (x, x') + (x)C(i) (x'), (57) 
V'(a;)V'(x')V'(a;") = C^^^ (a^, a;', a;") + C^^) (a;, xOC^^^ (x") + C^'^ (a;, a;")^^^) (x') 

(a:', a;")C(i) (x) + C^^' {x)C<-'^ ix')C^^^ {x"), (58) 



i.e., expression of a particular field average by the cumulants corresponds to its all possible factorisations. Character- 
isation of the field statistics by cumulants is more economical than characterisation by field averages. E.g., if "0(2;) is 
non-random, C^'^'^{x) = C^"^ = 0,n > 1; if ij{x) is random and Gaussian, its cumulants vanish for n > 2; 

whereas all field averages are non-zero even for non-random ip. 

There exists a general theorem (Mayer's first theorem) stating that the diagram expansion of the logarithm $(C) 
contains only connected diagrams, so that 

$(C) = exp[conn<i>(C)], (59) 

where conn$(C) is given by diagram series where all disconnected diagrams are dropped while connected ones retain 
their coefficients. Comparing this with the definition of the field cumulants, we see that 

conn$(C)=ln$(C)-^^C"C'("). (60) 

n=l ^' 

Hence (cf the examp les b elow) a diagram expansion for ;^C"C'''"'' contains all connected complete diagrams with 
exactly n C- vertices ( |52a| ), occurring with the coefficients they had in the series for $(C). Correspondingly, the 
expansion for C(") (xi,---,x„) contains all connected incomplete causal diagrams with exactly n free line outputs. 



C. Diagram structures corresponding to certain types of equations 

1. Emission of given sources 

To start with, we consider the simplest possible stochastic problem of radiation of a given random source, s{x) — 
sq{x). For simplicity (and also to stay within (pT])), we assume that it is Gaussian, and described by the cumulants 
so{x) — x^^'^H^) so{x)so{x') — so{x) so{x') — x^^''^^(x, a;') (conditioning on the full field is irrelevant for a given 
source). This problem is readily solved, 

^ix)=M^)+^'ix), (61) 
ip'{x) ^ I dx'G{x;x')so{x'), (62) 



(63a) 



C(2)(x,a;') = / dx"dx"'G{x;x")G{x'-x"')x^''''\x",x"') = < ^ , (63b) 



where ip'i^) is the emitted field. For the field cumulants we have. 



CW(x)=Vo(a:)+ J dx'G{x;x')x'^'''Hx') = [^} + l^ 




so that, 
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ln$(C)=CC« + ic^C(^) = {^} + {-^^} + i <r . (64) 



It is easy to see that these are indeed the three (and only three) connected complete diagrams that can be built of 
the available graphical elements -^Oj G, x^^''^^ and x'^'*^-*, with right coefhcients (cf the appendix). Hence expressions 
( |63| ) and (|6^) are exactly that one would find from the diagrammatic approach. Note, however, that the logic of the 
diagrammatic solution is reversed: relation ( |64|) is found summing all legitimate complete connected diagrams, while 
relations (|6^) follow. 

2. Linear susceptibility 

Consider now an equation with a non-zero linear susceptibility, x^^'^^ (x; x') ^ 0. For simplicity, we assume that 
the only other non-zero susceptibility is x^^''^''(^)i i-S-, we consider radiation of a given non-random source into a 
linear medium. The following considerations are nevertheless applicable in a general case of any set of non-zero 
susceptibilities. 

It is easy to see that all connected diagrams containing only the vertices x^^''^\x]x') and x^^'°-'(a:) are chains, 
1- 4(0 = { 3^ } + { ■ } + { ■ 3^ • } + ■ ■ + 

On their "past" end, the chains are terminated either by the line V'Oj 

{ ^ } - C^o, (66) 

{ ' } = CGx^'-'^V'o, (67) 

{ ' } ^ CGx^^'^^Gx^^'^Vo, (68) 

; (69) 

or by the vertex sq, 

{'^^7"}=CGx(^-"\ (70) 

{ • — — • } = cgx^^-^^gx^^-^), (71) 

{ • ' } = CGx(^'^)Gx(^-^)Gx(^-°\ (72) 

'■■ (73) 

The coefficients of the chains are all equal to one. (Formally, there is an additional class of connected diagrams, 
namely, loops, 



TrGx^'^'^ + ^TrGx^^'^^Gx'^^^) + iTrGx^^'^^Gx^^'^^Gx'^^^) + ■ • ■ • (74) 

They are zero due to causality conditions and regularisations. ) 

There are two possibilities as to how x^^'^^ can appear in a diagram: (i) between two Gs, and (ii) between a ipQ and 
G; it is the former that is responsible for infinite number of chains. Then, consider the sum of chains. 
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G"(x;x') = | l + l " ?| + |" — + (75) 

= G(x;x') + J dx"dx"'G{x-x")x^^'^\x"-x"')G{x"'-x') + ---. (76) 
It obeys an integral (Dyson) equation, 

G' = G + Gx^^'^'^G'. (77) 

Acting on it by the operator C, and using that CG = 1, we find, 

C'G' = 1, (78) 

where C = C—x^^''^^ ■ Thus partial summation of the chains in diagrams corresponds to shifting the linear susceptibility 
from the source to the "free" equation. 

Replacing G — * G' allows one to drop all diagrams containing a x^^''^'^ vertex between two G lines. As to the x*-^'^-* 
vertices placed between and Gs, these are only found in the combination, 

V'^,(x) = {—"} + { " } , (79) 

where the propagator is now G'. It is easy to see that CipQ — results in /C'V'o = 0. Hence by redefining the graphical 
notation, 

{— "}=^o(^), (80) 

= G'(x;x'), (81) 



= 0, (82) 

(and dropping the primes) we arrive at an equivalent problem where the linear susceptibility is included into the 
operator C. Below we always assume this to be the case. 

Note that instead of redefining the in-field, one could equally redefine the x^^'^"* vertex, 



(83) 

leaving unchanged. Then tpQ is no longer a solution to the "free" equation, but, firstly, when deriving diagrams 
using Eq. ( p8| ) this fact is irrelevant, secondly, £Vo = 0, and hence -CVo — 0, may anyway be not the case after 
regularisations. 



3. Nonlinear non- stochastic equations 

We now consider an equation with quadratic nonlinearity, x*-^'^-* ^ 0. We assume that this equation is non-stochastic, 
-^(2,0) _ x*-^'^-* = 0, and that x^^'^-* = (or included in the operator C). The available graphical elements hence are 
V'Oj G, x^^'°^ and x^^'^^ it is then easy to see that, with two exceptions, all connected diagrams that one can build 
using these elements are trees branching into the past: 

ln<i>(C) = {^} + {'^^7^} 
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The futuremost vertex in the trees is branching occurs at x^^''^^ vertices, and branches are ultimately terminated 
by tpo lines or x^^''^'' vertices. The series (jsj) contain all such trees, each with a coefScient 1/g, where g is the order 
of the group of symmetry of the tree. 



4-. Stochastic differential equations 

Consider firstly a special case when ;;^(™'") = for n > 1. Connected diagrams are then trees branching into the 
future: 



ln$ 




It is clear that diagrams with arbitrary number of ^ vertices may be found in this series. Hence, despite the fact that 
the local source conditioned on the full radiated field is Gaussian, the radiated field itself is non-Gaussian; moreover, 
it has non-zero cumulants C^"^ for all n. The number of connected diagrams is infinite, yet the diagram series retain 
certain "tameness": e.g., the number of diagrams contributing to each field cumulant, C*^"-*, is still finite. 

The final (and crucial) step leading to truly nontrivial series is combining nonlinearity and noise. Among the 
connected diagrams produced by relation (^), we find all trees, both "nonlinear" (Q) and "stochastic" (|8^); on the 
top of that, we find a totally new class of nctlikc diagrams which are often called diagrams with loops: 

ln$(C) = Trees+ i | <^^> « ]+U ^CD^ 




Now, the number of connected diagrams contributing to any field cumulant is infinite. 



IV. NONLINEAR STOCHASTIC RESPONSE PROBLEM 



The above CFT approach is a radiation problem, when one looks for properties of fields observable under given 
conditions; the quantum analog of it is, e.g., Glauber's quantum coherence theory [lO| ]. An alternative formulation 
is a response problem, when one is interested in the dependence of the field properties on external influences; in 
quantum mechanics, this is mainly associated with Kubo's linear response theory [ pT| . As was shown in [^,^, the 
formal diagram solution to the quantum nonlinear response problem has a structure identical to that of a classical 
nonlinear stochastic response problem. We therefore briefly outline how the above diagram expansions are adjusted 
to the "response" viewpoint. 

Following the spirit of Kubo's linear reaction theory we add a given external source to s{r,t) in the Eqs (|l|) 
and (^), s(r, t) — > s(r, t) + Scxt('', t), and consider how the properties of the field "ip depend on an infinitesimally small 
perturbation Scxt- We define nonlinear stochastic response functions as, 
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The response functions obey the causaHty condition, 

^(™,«) (2;^^ . . . ^a;^;a;'^^ . . . = 0, max(ti , • ■ •,*,„) < max(t'i, • • ■ ,t'^). 
For sufficiently small Sexti (and assuming we are not at a phase transition point) 



(87) 



00 ^ 

<i>(C,Sext)= 



m.n ) 



(89) 



m,n— 



where $(C,Scxt) is the characteristic functional of the field averages in presence of the external source, so that 

<i>(C) = 'i>(C,0). 

In diagram terms, adding the external source does not bring much new. One simply replaces x^^'*^' — > x^^''^'^ + Scxt: 



$(C,Soxt) = $(C)v(i,0)^y(l.0) 



(90) 



This means considering two types of x^^'"-* vertices and summing all diagrams thus emerging. If for simplicity we 
assume that x*-^'"-* = 0, the connection between the "radiation" and "response" series becomes trivial: To find the 
latter, take the former with 



= sext(a;). (91) 
Diagram expansions for the response functions are found by stripping the diagrams of the C and x''^'"' = Soxt vertices. 



R(^^^\x;x')^\ 



R^^'^-Xx^x'-x") 



R'^^'^\x;x',x") 




(92a) 



(92b) 



(92c) 



etc. It is easy to prove that the causality condition (88) holds not only for an _R(™'") as a whole, but also for each 
particular diagram in its expansion. 



V. CONSTRUCTING AN SDE FOR A GIVEN CAUSAL DIAGRAM SERIES 



Assume a diagram series is given (derived in a certain q- number approach, say). This series appears as a causal one, 
i.e., is generated by an expression like (48), and the causality conditions hold for all graphical elements. Formally, the 
generalised susceptibilities x^™'"^ then provide one with a complete and unambiguous description of the equivalent 
classical stochastic process. In practice, however, it would be more convenient to deal with an explicit SDE, written in 
terms of noise sources which are independent of the field. This leads us into what can be termed the inverse problem 
of the causal diagram techniques: How to explicitly write an SDE to which a causal diagram series corresponds? 

With no stochasticity present, the relation between a causal diagram series and the corresponding DE is straight- 
forward: A diagram series with only single-output vertices, x^^'"\ solves the equation. 
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°° 1 /■ 

£tp{x) = X! ^ / dxi ■ ■ ■ dxnX^^'^Kx; Xi, ■ ■ ■ , Xn)ip{xi) ■ ■ ■ 1p{Xn)- 



(93) 



Basically, any SDE should look identical to this equation, with the only difference that (some of) the susceptibilities 
are random. It is clear that, 



HO = $o(C), 



(94) 



where <I>(C) is the diagram series solving ( |93| ) regarded as an SDE, $o(C) is the series sol ving ( |93[ ) as a non-stochastic 
equation, and the upper bar here denotes averaging over the random susceptibilities in (|9^). The key to the inverse 
problem is in the relation between $o(C) and $(C)- 

Consider, to begin with, an example of a given random source. The SDE is then 



Cip{x) = so{x). 



(95) 

(96) 
(97) 



On the other hand, expanding <&o(C) in ^ diagram series and performing the averaging as required by (Q) in each of 
the diagrams separately, we find: 



For simplicity, we assume that sq is Gaussian with zero average, so{x) — and so{x)so{x') = x'^'^''^\x, x'). Then, 

$o(C)-exp{ • _ • }, 

$(C) = exp 




$(C) = 1 + { * 



1 + 0+- 




(98) 



We see that the graphical operation reflecting sq becoming random and Gaussian is a pairwise merging of the 
vertices into x^^'"-* vertices; diagrams with an odd number of the x'-^'"-' vertices become zero. Note that whereas itself 
the merging of the vertices reflects stochasticity, the fact that it is pairwise is clearly due to the Gaussian statistics: 
if So were non-Gaussian, then non-zero x*-™'"^ vertices would each result from merging of m x^'^'^'^ vertices. 
Consider now an SDE with a multiplicative noise, 



C^{x) 



dx'ri{x; x')ip{x') + s'(x). 



(99) 



where •q{x]x') is random and s'{x) contains non-stochastic terms (it may also contain other noise sources provided 
they are not correlated with rj). For simplicity, we again assume that r/(a;; x') is Gaussian and ry(x; x') = 0, so that it is 
specified by the average rj{xi] x[)ri{x2', a;^. The diagram series for equation ( |99| ) can be found by averaging those for 
an equation with a linear susceptibility x^^'^H^y — vi^y ^')- This is most intuitive if done directly in the generating 
expression (Ea): 



$(C) = exp 



/ s 



I exp {Ci)) S' (a|V') |a=o,i/-=V'o 

exp (CV') S' (a|V') la^o.v^'Ao' 



(100) 
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where S' {a\4>) contains vertices originating in s'(x), and the following graphical notation is used, 
{---"} =i,ix), 

{ 2ZZ } " 



X 

/ \ 



2 (2,2) ;2 



(101) 
(102) 
(103) 

(104) 



x'^'^'^\xi,X2;x[,x'2) = ?7(a;i;a;i)??(a;2;a;2) + '7(2^1; 4)'?(2;2; a;'i) 



(105) 



(The quantities a'"x''™'"'''0" ^Lre defined in Eq. (^3|). Note that we have departed from the case of only the five 
susceptibilities shown explicitly in Eqs. ( ^ ) being non-zero.) Hence the graphical representation of the linear 
susceptibility becoming random is pairs of the x^^'^' vertices merging into the x^^'^' vertices, e.g.. 




(106) 



Similarly, if a certain pair of susceptibilities (x''"'^'^^ and x^"'^'^^ say) become random and correlated (implying Gaussian), 
then averaging $0 we get, 



where 



m - exp [ -G- ) exp 



I 9 d 
exp — — G— — exp 



exp {(ip) S' {a\ip) |a=o,v=Vo 



^dtp da J ^12 
X exp (CV") S' {a\ip) \a=o,^=^o> 



48 



\ / 

X 

4 1/ \ 



1 



(107) 



(108) 
(109) 



(110) 
(111) 



X 

/ N 



2 (2.2) ;2 

a x^ ' V , 



(112) 



and 
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'°\xux2) = x(i'")(a;i)x(i'°)(x2), (113a) 
X^^''^' {xi,X2; x[,X2,x'^,x'4^) = Symmx^^'^) {xi;x[,x'2)x^'^''^') {x2;x'^,x'^), (113b) 
X^^^^\xuX2;x[,x'2) = Symmx(i-i)(xi;x;)x(i'i)(x2;4). (113c) 

Sy mm here denotes summation over all different terms obtained by permutations of the input and output arguments, 
cf (|l05| ). We see that whereas randomness results in self-mergings (x*^^'*^-* +x'^''''' x^^'^'' and +X*'^'^'' ~^ x'"^'"'^); 
correlations manifest themselves as cross-mergings (x^^'"' + X*^^'^'' ~^ x'^'^"*)- 

This way, if we take equation (p3|), derive a diagram series for it, and then assume that the susceptibilities x^^'"'' 
are random, the new series is found from the initial series by merging certain vertices in the diagrams. Hence in order 
to recover the initial non-stochastic equation, the merged vertices should be factorised into the products of the initial 



vertices. Formally, this means solving functional equations like (105) or (113), with given stochastic vertices, in order 
to find statistics of the susceptibilities in (|93|). This is nothing but the well-known problem of moments in probability 
theory. 

It is obvious that there exist causal diagram seri es that do n ot correspond to any SDE in the true meaning of the 
word. E.g., x^'^'^H^i x') — —S{x — x') would require so{x)so{x') — —5{x~x'); this is certainly impossible in probability 
theory. We therefore have to adopt the concept of pseudo-probability |lj] and consider pseudo-stochastic differential 
equations (PSDE) as well as stochastic ones. (This is anyway inevitable in quantum stochastics since the measure in 
Feynman path integrals is as a rule nonpositive.) Even with this generalisation, it is not clear if the inverse problem 
has a general solution consistent with the causality conditions for the susceptibilities. 

However, this clearly is the case for an important class of problems, namely, local Markovian problems. In terms of 
the susceptibilities, this means that they are non-zero only if all their arguments (both input and output) coincide. 
E.g., let x^^'^' (^ii ^2; a^i, 2:2) = X'^(2;i ~ x'i)S{x2 — X2)S{xi — X2), where x is a real constant, positive or negative. 
Then, the inverse problem is solved by x^^^'^H^', x') = ^/\x\1]ix)S{x — x'), where for x > ri{x) is a standardised 



Gaussian (5-correlated noise source, 'q{x)-q{x') = 5{x — x')\ whereas for x < 'q{x) is a standardised Gaussian 5- 
correlated pseudo- s tocha stic source, r]{x)ri{x') = —6{x — x'). Note that one can equally find an alternative solution 



using factorisation (112). This shows that a solution to an inverse problem is in general non- unique. 
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APPENDIX: 



We now show that relation (48) can formally be expanded in a series of the causal diagrams. To this end, we 
re- write it as. 



1 



plq\r\sW.u\ 



expl? 
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— 




I a—ip—0 7 



(Al) 



where 



and 



{ } = cv', (A3) 

|--<I| ^aV V (A4) 

^ (A5) 

The rest of the graphical notation used in ( |Al[ ) was introduced in Chapter 0. Note that we have returned to the case 
of only the five susceptibilities shown explicitly in E qs. ( |2l|) being non-zero. 



Rewriting the product of connected diagrams in (Al) as a disconnected diagram, the general term in the series 



(Al) is found to be, 

const X exp V {Dpqrstu} |a=v=o- (A6) 

Here {Dpqrstu} is a disconnected diagram with p + q + r + s + t + u connected components, where | * | is repeated 

p times, I *~ I is repeated q times, etc. To find a graphical representation of the operation exp V {Dpqrstu} \a=tii=o, 
we notice that, 

{-^pgrsiu} — dx\ • • • dXmdx^ • • • dx^ X 

i(;{xi) ■ ■ ■■4;{xm)a{x[) ■ • •a(a;^)S(a;i, ■ ■ ■ , x'^; xi, ■ ■ ■ ,Xm), (A7) 

where m = p + r + t + 2u, n = q + r + 2s + 2t + m, and S(x']^, • • • , x'^; xi, • • • , Xm) is a product of the generalised 
susceptibilies and ^s; a^i, • ■ • , Xm and x'^, • • ■ , are, respectively, input and output arguments, each belonging to a 
certain susceptibility factor (or to a C) in this product. S contains neither as nor ■(/;s, hence we have to understand 
only how exp V acts on the product of the as and ■i/'S. 

Consider firstly the operator expl?i, where 2?i = ^■00- It is readily seen that 

m 

Vi%l){xi) ■ ■ ■il){xm) = ^ V'(a;i) • ■ ■ i'{xi^i)^:o{xi)^:{xi+i) ■ ■ ■ijj{xrn)- (A8) 
1=1 

That is, Vi turns a product of m ips into m terms, in each of which one factor ip is replaced by tpo. Then, ^"Df 
will turn the same product into m(m — l)/2 terms, in each of which two factors V' are replaced by t/^qs; the factor 
^ exactly compensates for two possible orders of replacement leading to each term. It is now easy to realise (and 
prove by induction) that ^2?^ will turn the product of m t/js into different terms, in each of which k factors 

■0 are replaced by '0O- At fc = m, we shall find one term with all ips replaced by ipos, :^V^^p{xi) ■ ■ ■ 'ip{xm) = 
ipo{xi) ■ ■ ■ ipoixm), and V'^ipixi) ■ ■ ■ ip{xm) = 0,k> m. Thus, 

expPiV'(xi) • • • V(x„0 = ( 1 + 2?i + ^Vj + ■■■ + -^V'A V(xi) • • ■ V(a;„) 

\ z! m! J 

= Sum of all different products of k tpQS and m ~ k ipa, < k < m. (A9) 



Note that this result is based primarily on relation (AS); the rest of the derivation is mere combinatorics 



Similarly, how the operator expP2, where I?2 = '§^G-^, acts on a product of as and ips, is based primarily on the 
relation, 
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'D2ip{xi) ■ ■ ■ip{xm)a{x[) ■ • •a(a;^) = ^^G{xi; Xj)ip{xi) ■ • • -0(a;i-i)V'(a;;+i) • --ipixm) x 

1=1 j=i 

a{x[) ■ ■ ■ a(a;;-_i)a(x;.+i) ■ • • aix'J. (AlO) 

This relation reads: in the product ^{xi) ■ ■ ■ 'ip{xm)oi{x[) ■ ■ ■ a{x'^), select a pair 'ip{xi),a{Xj), and replace it by 
G{xi;x'j); 'D2ip{xi) ■ ■ ■ Tp{xm,)c({xi) ■ ■ ■ a{x'^) is then found as a sum of all dilferent terms obtained this way. Then, 
-g'D2'^{xi) ■ ■ ■ ip{x,n)a{x[) ■ ■ ■ a{x'^) is equal to the sum of all different terms, found by replacing k such pairs by the 
Gs in the initial product; if fc > min(m, n), one finds zero. Finally, 

exp'D2->p{xi) ■ ■ ■ 'ip{xm)a{x[) ■ ■ ■ a{x'„) 

— Sum of all different products of fc Gs, m — k and n — k as, < fc < min(m, n). (All) 

(This clearly follows the pattern of Wick's theorem for bosonic operators, replacing operator pairs by the "junction" 
G, cf m^.) Now there is no problem in understanding the action of the operator expP — exp(Pi + I?2): 

expVipixi) ■ ■ ■ ip{xm)a{x[) ■ ■ ■ a{x[^) 

= Sum of all different products of fc Gs, I ipoS, m — k — I ips and n — k as, 
0<k<n,0<l + k<m. 



(A12) 



The graphical recipe corresponding to this relation is quite trans paren t: (i) represent th e product 
tpixi) ■ ■ ■ ip{xm)oi{x'i) ■ ■ ■ a{x'„) as a disconnected diagram with m elements (101) and n elements ( |102| ), then (ii) 
write a sum of all different diagrams obtainable by replacing randomly chosen ■08 by ipoS and pairs a,')/' by Gs, 



{--:;}-{-■ 



(A13) 
(A14) 



For example. 



expP a{x)il;{x') = exp2? < ^ 



X X 



(A15) 



expV a{x)i/j{x')ip{x") — exp2? 




Then, setting a = = cancels all terms with dashed lines, e.g 

expP \ X t |a=i/'=0 



(A16) 



(Al7a) 



expP 



(A17b) 
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expV { x_ 

x' 



(Al7c) 



Thus, the graphical recipe for calculating eyi'p'D %p{xi) ■ ■ ■ ^j{xjn)a{x'i) ■ ■ ■ a{x'j^)\ci=,p=f) is that for calculating 
eyi'p'Dil){xi) ■ ■ ■ ^l){xm)oL{x'i) ■ ■ ■ a{x'^), plus (iii) retain only terms without dashed lines. Note that the number of 
resulting ip^s and Gs in all such terms is fixed {n and m — n, respectively), while their arguments vary from term to 
term. 



It is now obvious that to calculate (A6) graphically, one should remove pairs of dashed lines a^i}) in the diagram 
{Dpqrstu\i installing instead solid lines connecting inputs of the vertices to outputs, so that no free outputs remain; 
then, all remaining dashed ■(/^s should be replaced by solid V'oS, connected to inputs of the vertices. Term s w here the 
number of outputs in {Dpqrstu} exceeds the number of inputs do not contribute. Each terms in series (Al) is thus 
equal to either a sum of causal diagrams or zero, e.g., 



exp 2? I *" 

exp ^ I ~ 

iexpp{-_2 
expP| ^ I I I |a=v=o = ^expP I ^_ 



I |a=i/'=0 — 0, 
' I |a=0=O = I ' } 

" I |a=0=o = exp I? I 



|a=V'=0 




= 0=0 



(A18) 
(A19) 

(A20) 

(A21) 



(A22) 



The last exam ple shows that, in general, not all diagrams obtained from a particular {Dpqratu} are different (unlike 
in sums (A17)); this is due to possible symmetry of the vertex product S. All such equal diagrams count when 
calculating the coefhcient with which a diagram appears in the s eries . It is obvious that any given causal diagram, 
connected or disconnected, is produced by a unique term in series (Al), namely, by the one containing the same set of 
vertices as the diagram sought. It is then follows that all complete causal diagrams are found in the resulting diagram 
series. 
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